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ABSTRACT 

We study the non-thermal emission from old shell-type supernova remnants (SNRs) on 
the frame of a time-dependent model. In this model, the time-dependent non-thermal 
spectra of both primary electrons and protons as well as secondary electron/positron 
(e^) pairs can be calculated numerically by taking into account the evolution of the 
secondary pairs produced from proton-proton (p-p) interactions due to the acceler- 
ated protons collide with the ambient matter in an SNR. The multi- wavelength photon 
spectrum for a given SNR can be produced through leptonic processes such as elec- 
tron/positron synchrotron radiation, bremsstrahlung and inverse Compton scattering 
as well as hadronic interaction. Our results indicate that the non-thermal emission of 
the secondary pairs is becoming more and more prominent when the SNR ages in 
the radiative phase because the source of the primary electrons has been cut off and 
the electron synchrotron energy loss is significant for a radiative SNR, whereas the 
secondary pairs can be produced continuously for a long time in the phase due 
to the large energy loss time for the p-p interaction. We apply the model to two old 
SNRs, G8.7— 0.1 and G23.3— 0.3, and the predicted results can explain the observed 
multi-wavelength photon spectra for the two sources. 

Key words: radiation mechanisms: non-thermal - supernova remnants - ISM: indi- 
vidual: G8.7— 0.1, G23.3— 0.3 - gamma-rays: theory 



1 INTRODUCTION 

In recent years, various observational evidences indicate 
that particles in SNRs can be accelerated up to relativis- 
tic energies. For example. X-ray observations show that 
electrons can be accelerated to about 10 TeV in the shell 
of an SNR (Koyama et al. 1995; Reynolds 1996; Bamba 
et al. 2003, 2005a, 2005b). TeV 7-ray observations indi- 
cate that both electrons and protons can be accelerated 
to several tens TeV in young shell-type SNRs such as 
RX J1713. 7-3946 (Enomoto et al. 2002; Aharonian et al. 
2006a) and RX J0852.0-4622 (Katagiri et al. 2005; Aharo- 
nian et al. 2005). More recently, a survey of the inner part 
of the Galactic Plane in VHE 7-ray band has been per- 
formed with HESS Cherenkov telescope system. Aharonian 
et al. (2006b) presented detailed spectral and morphologi- 
cal information of some new sources discovered in the HESS 
survey along with some discussions on possible counterparts 
in other wavelength bands. Some sources discovered in the 
survey can be associated with old SNRs. For example, HESS 
J1804-216 has a flux of nearly 25% of the flux of the Crab 
Nebula with a photon index of about 2.72. The source can 
be Eissociated with either the shell-type SNR G8.7— 2.1 or 
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the young Vela-like pulsar PSR J1803— 2137 (Aharonian et 
al. 2006b). Furthermore, Fatuzzo, Meha, & Crocker (2006) 
found that the SNR G8.7— 2.1 can plausibly account for all 
known radiative characteristics of HESS J1804— 216. An- 
other case is HESS J1834-087, which is an extended TeV 
source discovered in the HESS survey of the inner Galaxy in 
VHE 7-ray band. The source, spatially coincident with the 
SNR G23.3-0.3, has a 7-ray flux of 8% of the flux of the 
Crab Nebula above 200 GeV, and the 7-ray spectrum can be 
expressed by a power law with a photon index of 2.45 ±0.16 
(Aharonian et al. 2006b). The new XMM-Newton observa- 
tion reveals diffuse X-ray emission within the HESS source 
and suggests an association between the diffuse X-rays and 
the VHE 7-rays, moreover, G23.3— 0.3 is estimated to be an 
old SNR with a distance of 4 ± 0.2 kpc and an age of ~ 10^ 
yr (Tian et al. 2007). 

Theoretically, it is believed that the non-thermal pho- 
ton spectrum for a given SNR is produced through elec- 
tron/positron synchrotron radiation, bremsstrahlung and 
inverse Compton scattering as well as tt" decay from the 
p-p interactions. The question is, however, whether the ori- 
gin of the non-thermal photons for young SNRs is different 
from that for old ones. Sturner et al. (1997) presented a 
temporal evolution model to calculate the non-thermal par- 
ticle and photon spectra from shell-type SNRs, and Zhang & 
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Fang (2007) modified the model to investigate the possible 
hadronic contributions to TeV 7-ray emission from young 
SNRs. For the young SNRs, the secondary particles pro- 
duced via the p-p interaction have a negligible contribution 
to the non-thermal photon spectrum (e.g. Zhang & Fang 
2007). On the other hand, Yamazaki et al. (2006) studied 
the emission from old SNRs and found that proton accelera- 
tion could be efficient enough to emit TeV 7-rays both at the 
shock of the SNR and at the giant molecular cloud. They 



presented that the energy flux ratio R^ 



ToV/X 



can be more 



than lO'^ for some SNRs, and such sources may be the origins 
of the recently discovered unidentified TeV sources. For the 
old SNRs, the energies of the primary electrons accelerated 
by the shock wave are limited by the synchrotron cooling, 
then the roll-off energy of the synchrotron emission of the 
primary electrons is much smaller than that for the young 
SNRs (Sturner et al. 1997; Yamazaki et al. 2006; Zhang & 
Fang 2007). However, the secondary pairs can be pro- 
duced from the p-p collisions in an SNR for a long time due 
to the large energy loss time for the p-p interaction. There- 
fore, the contributions of secondary electrons and positrons 
to the non-thermal emission can become prominent for the 
old SNRs. 

Motivated by above discussions, we study the non- 
thermal photon spectra from old SNRs based on the time- 
dependent model in which the contribution of the secondary 
pairs is included (see §2.2). It should be noted that our 
model is different from that given by Yamazaki et al. (2006) , 
the latter assumes a steady particle distribution. This paper 
is organized as follows. The temporally evolving model of 
the non-thermal particle and photon spectra are discussed 
in Section (2] including the SNR evolution, shock accelera- 
tion, and various photon production mechanisms involved in 
the model. The results from the applications of the model 
to two old SNRs, G8. 7-0.1 and G23. 3-0.3, are shown in 
Section (3] The main conclusions and some discussions are 
given in Section [l] 



TEMPORAL MODEL FOR 

MULTI- WAVELENGTH NON-THERMAL 

EMISSION FROM SNRS 



ISturner et al. I (|l997l l presented a temporally-evolving 
model for the non-thermal particle and photon spectra at 
different stages in the lifetime of a shell-type SNR. Zhang 
& Fang (2007) modified the model to explain the multi- 
wavelength observations of the young SNRs and shown 
that the TeV 7-rays from the two young shell-type SNRs 
RX J1713.7-3946 and RX J0852. 0-4622 have hadronic ori- 
gin. In this paper, we include the evolution of secondary 
pairs to model the multi-wavelength non-thermal emission 
from old SNRs since the secondary pairs can have sig- 
nificant contributions to the non-thermal emission from old 
SNRs. 



2.1 SNR Evolution and Shock Acceleration 

In the analytical model of the shock dynamics of an SNR ex- 
panding into a uniform ambient medium with density no , the 
SNR evolves through three stages: the free expansion stage, 
the Sedov stage, and the radiative stage (Sturner et al. 1997; 



Yamazaki et al. 2006) . Assuming the initial explosion energy 
of the SNR \a E = lO^^Ssi ergs, and the initial velocity of 
the shock is vq, the free expansion stage ends when t = 
tscd = {3E/2nmHnoviy^^ « 2.1 x lO'' {Esi/noy^^'vo^^^ yr. 



and the Sedov stage ends at t = t,. 



4.0 X 



4/17 -9/17 



yr (Blondin et al. 1998), where no = fmisM, nisM is the hy- 
drogen density in the local interstellar medium, fi = 1.4 is 
the mean atomic weight of the interstellar medium assuming 
1 helium atom for every 10 hydrogen atoms, and mn is the 
mass of hydrogen. The shock velocity Vs{t) can be expressed 
as (e.g. Yamazaki et al. 2006) 



vs{t) 



Vo 
Vo 



*Sod/ 



-3/5 
-3/5 



-2/3 



t < tScd 

^Scd ^ t < ti-ad (1) 
trad ^ t . 



The evolution of the shock radius is estimated by Rs{t) = 
/ Vs{t)dt and we assume vo = 10^ cm s~^ in this paper. It 
should be noted that the accumulated energy in non-thermal 
protons for an SNR is usually smaller than 10% of the initial 
explosion energy (see details in §2.1), and the energy dimin- 
ishes very slowly in the radiative phase when the thermal 
energy in the SNR experiences significant loss. At late times, 
the relativistic protons begin to contribute significantly to 
the post-shock pressure when most thermal energy of the 
SNR is drained and then the expansion law should be re- 
vised correspondingly. However, the expansion law will be 
infiuenced significantly by the accelerated protons at very 
late times and the non-thermal emission from old SNRs does 
not depend heavily on the expansion law. Therefore, we still 
assume the expansion law expressed by Eq. ((T} is always 
valid in this paper. 

It is generally believed that the non-thermal charged 
particles are produced through diffusive shock acceleration 
(Blandford & Eichler 1987; Elhson, Jones & Reynolds 1990). 
Neglecting the nonlinear effects and assuming a compression 
ratio near 4, the particle spectrum from the diffusive accel- 
eration process is a power law with spectral index a ~ 2. 
For the particles accelerated in the shock of an SNR, the 
power law particle spectra do not extend to infinite energy 
and must be truncated by three mechanisms, namely, the 
finite age of the SNR, energy-loss processes, and free escape 
from the shock region (Sturner et al. 1997). The maximum 
kinetic energy of electrons and protons, E^^ max and -Ep, max, 
can be calculated by using Eqs. (3), (4), (5) in Sturner et 
al. (1997), and we neglect the exact formulae here. Zhang 
& Fang (2007) studied the non-thermal emission from three 
young SNRs, and concluded that the model in the paper 
can explain the multi-wavelength observations for the three 
young SNRs with the maximum wavelength of MHD turbu- 
lence Amax = 2 X lO'^^'cm and other appropriate parameters. 
Amax is also set to 2 x lO^^cm in this paper. 

It should be pointed out that a time-dependent ki- 
netic equation could be solved in order to obtain the en- 
ergy spectrum of the accelerated particles (Berezhko et 
al. 1996; Sturner at al. 1997). Here, following Sturner at 
al. (1997), we approximate the volume-averaged emissivity, 
Q{E, t) = dN/dVdtdE, of the shock accelerated electrons 
and protons by 



Qr{E,t) = Q'iG{t)[E{E + 2m^d')] 



-[(a + l)/2] 
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x{E + rUeC?) exp(- 



E 



At)' 



(2) 



and 

QfiE,t) 



2 E 
x{E + mpC )exp{-- — ) , 

-C^p, max [J' ) 



(3) 



respectively, where E is the particle kinetic energy, and nie 
and nip are electron mass and proton mass, the index a ~ 
2.0, and we use a = 2.0 in this paper. The function G{t) is 
given by (Sturner at al. 1997) 



G{t) 



[-RsNR(iScd)/-RsNa(t)] t ^ trad 
t > trad 



(4) 



Factors and Qp are used to normalize the particle spec- 
tra so that the total amount of the kinetic energy contained 
in the injected electrons and protons is i?par = rjMejVQ/2, 
where 77 ~ 0.1 represents the efficiency that the kinetic en- 
ergy of the ejecta with initial mass Afcj and initial velocity 
Ho is converted into the kinetic energy of both the electrons 
and the protons, 



So 



dtVsMt){ f " dEEQT{E,t) 
Jo 



dEEQf{E,t)) 



(5) 



and VsNR(f) = 47rJ?iNj,(t)/3 is the SNR volume. In Eq. ®, 
we need to introduce the parameter Kep = Q^/Qp in order 
to determine and Qp. Obviously, these two quantities de- 
pend on rj and Kcp. The ratio of electrons to protons, Kcp, is 
a key parameter when modeling the multi-wavelength non- 
thermal emission from an SNR. The measured ratio from 
the cosmic ray observations at the Earth is about 0.013 at 
10 GeV (Gassier 1990). However, the measured value can 
not represent the ratio for a given SNR due to the cosmic 
ray transport effect in the Galaxy. Sturner et al. (1997) ob- 
tained Kcp ~ 0.6 by assuming that the kinetic energy of the 
injected electrons is the same as that of the injected pro- 
tons. Yamazaki et al. (2006) studied the TeV emission from 
the old SNRs in the steady state, they used the observed 
ratio of TeV 7-ray (1-10 TeV) to X-ray (2-10 keV) energy 



flux, R 



, -ftTeV/X, 



for some SNRs to test the values of Kcp, and 



2.2 Temporal Evolution of Particle Energy 
Distributions 

Following Sturner et al. (1997), assuming that an SNR inte- 
rior is homogeneous, with a constant density tisnr = 4nisM 
and a magnetic field strength Bsnr = 4_Bism, ne{Ee,t) 
and np{Ee,t) are used to represent the differential densities 
of the accelerated electrons and protons, respectively. We 
can use Qe and Qp to calculate the direction- and volume- 
averaged electron intensity Js{Ee,t) = (c/3/47r)ne(-Be,t) and 
proton intensity Jp{Ep,t) = {cl3/4:n)np{Ep,t) at each mo- 
ment during the SNR lifetime. Fokker-Planck equations for 
both electrons and protons in energy space are used to find 
solutions, which are given by (Sturner et al. 1997) 



dne{Ee,t) 
di 



-^^ETn4E.,t)] 

+ l-^[D{E,,t)n4E,,t)] 



+ Qe{E,,t)- 



ne{Ee,t) 



(6) 



and 

dnp{Ep,t) 

Ft 



+ ^^[D{Ep,t)np{Ep,t)] 



+ Qp{Ep,t) 



np(Ep,t) 



(7) 



respectively, where the terms on the right-hand sides in 
Eqs. and (0 represent systematic energy losses, diffu- 
sion in energy space, the particle source function and catas- 
trophic energy loss. Following Sturner et al. (1997), we solve 
the above equations using a Crank-Nicholson finite differ- 
ence scheme. In Eqs. (|6} and ©, El°\ D{Ee,t), El°^ and 
D{Ep,t) can be calculated with the formulae in §2.2 in 
Zhang & Fang (2007), the exact expressions are neglected 
in this paper. 

The secondary electrons and positrons can be produced 
when the accelerated protons collide with the ambient mat- 
ter via the p-p interaction. Those secondary particles evolve 
with the SNR, and can contribute to the multi-waveband 
non-thermal emission from the SNR as the primary elec- 
trons. So the source term for electron can be represented 
as 



found that the predicted i?ToV/x for young SNRs is similar q^^^^-^-j ^ gpri _^ gscc (^^ _^ qbcc ^-^^ 



to the observed one when Kep = 1 x 10~ , and that the pre- 
dicted values of -RtoV/x for the old SNRs are not changed 
even if Kep varies more than one order of magnitude because 
the primary electrons do not contribute to TeV and X-ray 
emission in their model. Moreover, Berezhko, Ksenofontov & 
Volk (2006) investigated the properties of Kepler's SNR us- 
ing the nonlinear kinetic theory of cosmic ray acceleration in 
the SNR. They treated Kep as a parameter quantitatively 
determined by comparison with the synchrotron observa- 



where 



Q 



{E,t) = 471 fippUsNR J dEpJp(Ep,t)'^'^^^' 



, Er, 



dE, 



QT-{E,t) =47r^ippnsNR / dEpJp{Ep,t) 



^ da(E^- , Ep 



dE, 



(8) 



(9) 



(10) 



tions and obtained that A'c 



10"" for the SNR. In fact. 



the value of Kep for a given source is usually uncertain now, 
and it can be limited by comparison of the resulting spec- 
trum with the multi-wavelength observations for an SNR 
(Zhang & Fang 2007). In this paper, 77 is set to 0.1 and 
Kep is treated as a parameter limited by comparison of the 
calculated results with the observations for a given SNR. 



Mpp = 1-45 is a enhancement factor for collisions involving 
heavy nuclei in an SNR (Dermer 1996; Sturner et al. 1997), 
da{Ee-,Ep)/dEe- and da{Ee+ , Ep)/dE^+ are the differen- 
tial cross section for electrons and positrons produced via 
p-p interaction, respectively. Since the accelerated protons 
colliding with the ambient medium (i.e. p-p interaction) 
can produce pions (tt^ and n"), in which a proton loses 
~ 1/3 of its energy per pion-producing collision, there exists 
the catastrophic loss for the protons. Sturner et al. (1997) 
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treated this process as an escape from the system, and the 
timescale for this catastrophic loss is 



rpion(-Bp) = [c/3pnsNRO-pp] 



(11) 



where cTpp is the inelastic cross section for p-p interaction, 
which is (Kelner et al. 2006) 



crpp(£;p) = (34.3 + 1.88L + 0.25L^ 



1 - 



mb, (12) 



where L = ln(ii'p/lTeV), -Bth = nip + 2m,r + m^/2mp is the 
threshold energy of production of 7r-mesons. 



2.3 Photon Production 

We can solve the Eqs. © and (0 to get the volume- and 
direction-averaged intensities for electrons and protons at a 
given time, and then the multi-waveband non-thermal pho- 
ton spectra can be calculated. The non-thermal radiation 
processes of the accelerated particles involved in an SNR 
are synchrotron radiation, bremsstrahlung, inverse Comp- 
ton scattering for leptons including electrons and positrons, 
and the p-p interaction for protons. 

Accelerated electrons can emit photons via synchrotron 
radiation as interacting with the ambient magnetic field in 
an SNR. The emissivity for synchrotron radiation is given 
by 

Qsyn{E-,,t) = — — — / d9 sm 6 

\ hE^nieC-^ ) Jo 

'dSeJe(£e,t)i^('^V (13) 



J E, 



where Q is the electron pitch angle, Ec = 4.2 x 

„2 , 



10^ ft-Bs Affile sin ^, and 
F{y)=y dzKsrAz) 

J V 



(14) 



with y = E-y/Ec, where if 5/3 is a modified Bessel function 
of order 5/3. 

For the electron-electron and electron-nucleon 
bremsstrahlung radiation, the emissivity is (Sturner et 
al. 1997) 



]hi-cm{E^,t) = 47rnsNRA^'" 

w , ^ da 

X / dEeJe{Ee,t) 



dE. 



(Ifi) 



where A^e and Ag;, are correction factors for the presence 
of helium (Age = 1-2 and Agg = 1.4), drje-e/dE^ and 
dae-p/dEj represent the differential electron-electron and 
electron-nucleon bremsstrahlung cross sections. Here, we use 
the approximate formulae given in Baring et al. (1999). 

We use the full Klein-Nishina cross section for the rela- 
tivistic electrons to tackle the inverse Compton process, the 
emissivity can be expressed as 



roc rEerna^x 

Qcomp, j(-B^,t) = 47r / denj{e,r) / dEe 
xJ4Ee,t)F{e,E^,Ee) , 



(16) 



where, e is the target photon energy, nj is the number den- 
sity of soft photon component j with energy density Uj and 
temperature Tj and is given by 



15C/i 



1 



(17) 



(^fcTj)4exp(e/fcTj) 

-Ee.thrcsh = [E-, + (E?, + E^{meC^f /ey^^]/2 is the lowest en- 
ergy that electrons can scatter a target photon with energy 
e to energy E^. Function F{e, E-y, Ee) is given by 

3aT 1 



F(e,E^,E, 



AiE.Jmc'^Y e 
2glng+(l + 2g)(l-g) 



(rg)'(i- 



(18) 



2(1 + rg) 

with r = 4:e{Ee/mc^)/mc^, q = Ei/r{l - Ei) with Ei = 
E^/Ee and l/A{Ee/mc^) < q < 1. 

Kamae et al. (2006) presented a series of accurate, 
convenient parameterized formulae to calculate the spec- 
tra of stable secondary particles (7, e*, Ue, Ve, u^, Vp) 
produced in the p-p interactions, which greatly facilitate 
calculations involving the p-p interaction in astronomical 
environments. The formulae were derived from the up-to- 
date p-p interaction model given by Kamae et al. (2005), 
which incorporates the logarithmically rising inelastic cross 
section, the diffraction dissociation process, and the Feyn- 
man scaling violation. The functional formula reproducing 
the secondary particle spectra include the non-diffractive, 
diffractive, and resonance-excitation components. For the 
nondiffractive process, the differential inclusive cross section 
to produce a secondary particle is given as 

S< 



= Fm{x)Fm,iii{x) 



(19) 



dln(£;scc) 

where E^cc is the energy of the secondary particle, x — 
Escc/GeV, Ftiu{x) is the formula representing the non- 
diffractive cross section, and -F'ND,ki(a:) is the formula to 
approximately enforce the energy-momentum conservation 
limits: 



Fnd{x) = aoexp{—ai{x — aa + a2{x — as)^)'^) 
+ 04 exp(— a5(a- ~ as + a^ix — as)^ 
+ a7{x - ag)^)^), 



-Fnd,] 



exp [WND,l(imin - x)] + 1 
1 



(20) 



(21) 



exp [iyND,h(3:^ - i/max)] + 1 

where Lmin and Lmax are the lower and upper kinematic 
limits imposed and Wnd,i and Wnd.h are the widths of the 
kinematic cutoffs, Lmin ~ —2.6 for all secondary particles, 
and the other parameters are given in Table [1] (Kamae et al. 
2006). 

The differential inclusive cross section to produce a sec- 
ondary particle for the diffractive process can be represented 

as 

rfo"diff(-Bscc) 



FdiH(x)FM{x) , 



(22) 



d\n{Eaec) 

where Fdia{x) represents the diffractive cross section, and 
-^111(2;) is a function to enforce the energy-momentum con- 
servation: 

i^diff(a;) = boexp{-bi{{x-b2)/{l + ba{x-b2))f) + 

64 expi-b^iix - 66)/(l + b7{x- be))f), (23) 
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Figure 1. The isotropic intensities of the primary electrons (dotted Unc), the primary protons (solid line), and the secondary pairs 
(dashed line) at (a) 2000, (b) 8000, (c) 15000, and (d) 30000 yr, respectively, for A4j = l.SM©, nisM = 10 cm'^, Bism = 5nG, 
= 0.001. 



Table 1. Kinematic limit parameters for the nondiffractive pro- 
cess (Kamae et al. 2006). 





^max 


l^ND.l 


H^ND,h 


7 


0.961ogio(Tp) 


15 


44 


e~ 


0.961ogio(Tp) 


20 


45 


e+ 


0.941ogio(Tp) 


15 


47 


l^e 


0.981ogio(Tp) 


15 


42 


Pe 


0.981ogio(Tp) 


15 


40 




0.941ogio(Tp) 


20 


45 


Pm 


0.981ogio(Tp) 


15 


40 


1 



exp (Wdiff(a:: - imax)) + 1 ' 



(24) 



here Wmtt = 75, L„,ax = logio(rp/GeV). 

For the resonance-excitation processes, the function is 



(25) 



darasiEs, 

dln(£sc. 

Frcs{x) represents the cross section, and Fu{x), which is 
same as that for the diffraction process, enforces the energy- 
momentum conservation: 



Fresix) = Coexp{-Cl{{x - C2)/{1 + C3{x - C2) 
+ C4{x - C2)^))^). 



(26) 



A renormalization factor, r{Tp), must be multipUed to 
the final spectrum to ensure that the parameterized model 
reproduces the experimental tt" multiplicity after the read- 
justment in the resonance-excitation region of Tp. The ex- 



act expressions for r(Tp) and other parameters involved in 
the formulae are shown in Kamae et al. (2006). Finally, the 
differential cross section for the secondary particles can be 
expressed as 



da[Eaac, Ep 
dEsec 



da-m{Escc, Ep) rfcdiff (-Escc, Ep 



+ 



+ 



dEscc dEsec 
rfcrA(1232)(-E'sec, Ep) 



dE,, 



+ 



da 



ros(1600) 



{Escc, Ep) 



dE„ 



(27) 



The resulting particle intensities calculated with the 
model are shown in Fig. [1] The parameters we choose are 
Moj = 1.5Mo, nisM = 10cm"^ Bism = 5/iG, K^p = 0.001, 
for ages of 2000, 8000, 15000, 30000 yr. For these parameters, 
the primary particle sources turn off at f = frad = 1.01 x lO'* 
yr. The turnover at the low energy region for both electrons 
and protons is due to Coulomb loss, and the high-energy 
turnover for electrons is due to synchrotron loss. More im- 
portantly, the primary electrons dominate the secondary 
pairs for a young age, however, the secondary pairs be- 
come more and more important and eventually surpass the 
primary electrons as the age of the SNR increases, the rea- 
son is that the sources of the primary particles are cut off 
in the radiative phase of the SNR and the primary electrons 
experience strong synchrotron energy loss, however, the en- 
ergy loss for the protons is negligible due to the small cross 
section for the p-p interaction, and then the secondary par- 
ticles can be produced continuously for a long time. 

We show the resulting multi-wavelength non-thermal 
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Figure 2. The resulting calculated multi-wavelengtii spectra at (a) 2000, (b) 8000, (c) 15000, and (d) 30000 yr, respectively. The distance 
is set to 2 kpc and the other parameters are same as those in Fig. [T] The emission due to synchrotron (dashed line) , bremsstrahlung 
(dotted line), Compton scattering (dash-dotted line) of the primary electrons (thin) and the secondary pairs (thick), and the p-p 
interaction (solid line) are shown in each panel. 



spectra in Fig. [2] The distance to the SNR is chosen to be 
2 kpc, and other parameters are same as those in Fig. [1] 
It is obvious that the emission from radio to X-ray band 
for the SNR comes from the synchrotron radiation of lep- 
tons, and the high-energy photons can from bremsstrahlung 
and inverse Compton scattering of leptons, and meson de- 
cay due to the hadronic interactions. More importantly, (1) 
the roll-off energy of the synchrotron radiation of the pri- 
mary electrons deceases quickly when the SNR ages in the 
radiative phase, however, the synchrotron emission of the 
secondary pairs become more and more important and 
eventually surpass the synchrotron radiation from the pri- 
mary electrons when the SNR ages due to the secondary 
pairs can be produced for a long time via the p-p interaction; 
(2) most VHE 7-rays from the SNR are of hadronic origin 
with those parameters, and the hadronic emission from an 
SNR is similar during a long time when the SNR is in the ra- 
diative phase because the energy loss rate for the contained 
protons is small. Of course, the photons with energy smaller 
than 10 MeV from some old SNRs mainly come from the 
radiation of the secondary pairs. 



3 APPLICATION TO OLD SNRS 
3.1 SNR G8.7 0.1 

Aharonian et al. (2006b) reported a survey of the inner part 
of the Galactic Plane in very high energy 7-rays observed 
with HESS Cherekov telescope system, and fourteen previ- 
ously unknown sources were detected at a significance level 
greater than 4cr after accounting for all trials involved in the 
search. A source detected in the survey is HESS J1804— 216, 
which is the brightest one of the new sources, with a steep 
photon index of 2.72 ± 0.12 and a flux of nearly 25% of the 
flux from the Crab Nebula above 200 GeV. The source can 
be associated with the south-western part of the shell of 
the SNR G8. 7-0.1 of radius 26 arc minutes with flux 80 Jy 
at 1 GHz (Green 2006), or the young Vela-like pulsar PSR 
J1803— 2137 with a spin-down age of 16000 yr. Fatuzzo et al. 
(2006) presented a simple model for the TeV emission from 
HESS J1804-216, they found that the SNR G8.7-0.1 can 
plausibly account for all the known radiative characteristics 
of HESS J1804-216, and therefore the SNR G8.7-0.1 was 
probably the source of the TeV photons originating from the 
direction of HESS J1804-216. 
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Figure 3. Comparisons of the model results with the multi-band observations for the SNR G8.7— 0.1. The radiation from the primary 
electrons (dashed line), the primary protons (solid line), and the secondary pairs (dotted line) with data points from Odegard (1986), 
Kassim & Weiler (1990), Green (2006) (radio), Cui & Konopelko (2006) (X-ray), Aharonian et al. (2006b) (HESS) are indicated in each 
panel. The parameters are a distance of 6kpc, an age of 15000 yr, M^j = 1.5 M0, njsM = 8cm~'^, -Bism = 6/xG, (a) Kcp = 0.001, (b) 
K^p = 0.01, (c) Kep = 0.1, (d) K^p = 0.2. 



By associating G8.7— 0.1 with coincident HII regions, 
Kassim & Weiler (1990) estimated a distance to the SNR of 
6 ± 1 kpc, with angular size of ~ 50 arcmin, a physical size 
of ~ 80pc, and an age of 15000 yr. Cui & Konopelko (2006) 
reported the high-resolution X-ray observations taken with 
the Chandra X-Ray Observatory of the field that contains 
the TeV 7-ray source HESS J1804-216, and a total of 11 
discrete sources were detected. Among those sources, only 
one, CXOU J180351.4-213707, is the most probable X-ray 
counterpart of HESS J1804— 216, is significantly extended 
with photon index 1.21q4. 

We now assume that the sources CXOU J180351.4 
- 213707 and HESS J1804-216 associate with the SNR 
G8.7— 0.1. To make the resulting spectrum consistent with 
the observations by Chandra with a low energy fiux and a 
photon index about 1.2, the SNR must be in the radiative 
phase when the roll-off energy of the synchrotron radiation 
of the primary electrons decreases quickly and the X-rays 
detected by Chandra X-Ray Observatory are most proba- 
bly from the bremsstrahlung radiation or inverse Compton 
scattering of leptons, then nisM should be greater than ~ 5 
cm~"^ for an age of 15000 yr. 

We now describe the detailed process in modeling the 
multi-waveband spectrum for the SNR. First of all, we inves- 
tigate the influence of K^p on the final spectra with param- 
eters a distance of 6 kpc, an age of 15000 yr, Mcj = 1.5 Mq, 
Bism = 6/iG, and msM = 8cm~=' for the SNR G8. 7-0.1. 
The calculated results are shown in Fig. [S] The spectrum 



in the X-ray band depends heavily on K^p, and the value 
of Kep must be around 0.1 in order to match the observed 
data in X-ray band for the SNR. 

Fig.|3]shows the resulting spectra from the SNR for dif- 
ferent nisM and ages. For an age of 15000 yr, the roll-off 
energy of the synchrotron emission is still high for nisM = 
6cm~'^, which conflicts with the Chandra observation, and 
the level of the TeV emission is greater than the observed one 
by HESS both for tiism = 6cm~'^ and for nisM ~ 8cm~'^ 
(see penal (a) and penal (b) in Fig. 3)). The model results 
can fit the multi-band observations well with a lager age of 
18000 yr, moreover, the calculated spectrum in the band 0.1 
- 1 TeV for nisM = 6cm~^ fit the TeV observations better 
than that for tiism = 8cm~^ (see penal (c) and penal (d) in 
Fig.|4ll. Furthermore, the contributions to the final spectrum 
of different radiation components with the same parameters 
as those in penal (c) of Fig. |4] are shown in Fig[5] Obvi- 
ously from FigO for the SNR, (1) the emission from the 
primary electrons dominates that from the secondary 
pairs in the entire energy band except in the narrow soft X- 
ray band around 0.5 keV; (2) the detected radio emission is 
mainly from the synchrotron radiation of the primary elec- 
trons whereas the X-rays observed with Chandra are primar- 
ily produced via the bremsstrahlung of these electrons; (3) 
the TeV photons with energies < 1 TeV are primarily from 
both the bremsstrahlung of the primary electrons and the 
p-p interaction of the primary protons, however, those with 
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Figure 4. Comparisons of the model spectra witli tlie observations for tlie SNR G8.7— 0.1 with K^p = 0.15, (a) nisM = 6cm~^, an age 
of 15000 yr; (b) nigM = 8cm~'^, an age of 15000 yr; (c) nisM = Gcm~^, an age of 18000 yr; (b) nisM = 8cm~^, an age of 18000 yr. The 
whole emission from the primary particles and the secondary pairs are indicated by thin solid line and thick solid line respectively. 
Others are same as those in Fig. \3\ 
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Figure 5. The comparison of detailed spectra calculated by the 
model with the observations for the SNR G8.7— 0.1. The parame- 
ters are same as those for panel (c) in Fig.|4] Synchrotron (dashed 
line), bremsstrahlung (dotted line), Compton (dash-dotted line) 
of the primary electrons (thin) and the secondary pairs (thick), 
meson decay due to the p-p interaction (black solid line) are in- 
dicated in the figure. 



higher energies are mainly from the p-p interactions due to 
the high-energy protons collide with the ambient matter. 



3.2 SNR G23.3-0.3 

HESS J1834-087 is another extended TeV source recently 
discovered in the HESS survey of the inner Galaxy in VHE 
7-rays. The source, which is spatially coincident with the 
SNR 023. 3-0.3 (W41), has a size of 12 arc minutes and a 
7-ray flux of 8% of the flux from the Crab Nebula above 
200 GeV, with a photon index of 2.45 ± 0.16 (Aharonian et 
al. 2006b). Moreover, Albert et al. (2006) presented the ob- 
servations of HESS J1834-087 with the Major Atmospheric 
Gamma Imaging Cherenkov telescope {MAGIC), and the 
differential 7-ray flux can be expressed as dN / dAdtdE = 
(3.7 ±0.6) X 10~^^(£;/TeV)-^ '^±°-^ cm'^ s"^ TeV^ 

G23.3— 0.3 is an asymmetric shell-type SNR with radio 
and X-ray photons having been detected. In radio band, the 
SNR has a spectral index of 0.5 and a flux of 70 Jy at 1 
GHz (Green 2006). Moreover, Tian et al. (2007) reported 
new H I observations from the VLA Galactic Plane System 
and a new XMM-Newton observation for HESS J1834-087. 
They concluded that G23.3— 0.3 is an old SNR with a dis- 
tance of 4 ± 0.2 kpc and an age of ~ 10^ yr. Furthermore, a 
density of 6 cm~^ was estimated around the SNR from the 
VGPs H I line emission associated with the SNR (Tian et al. 
2007). The new XMM-Newton observation indicated diffuse 
X-ray emission within the HESS source and suggested an 
association between the X-ray and the VHE 7-ray emission. 
The intrinsic X-ray spectrum from 2 to 10 keV could be 
characterized as a heavily absorbed power law with a flux 
of 7.0 X 10"^^ erg s"^ cm"^, and a photon index of 2.01q g 
(see details in Tian et al. 2007). 



Non-thermal emission from old supernova remnants 9 




10 

10" 



10 " 10 " 10"^ 10" 10 ' 10"" 10' 10' 
Photon Energy (MeV) 



10 " 10" 10' 10" 10' 10' 10' 10" 
Photon Energy (MeV) 




10"'" 10"'" 10"' 10" 10"' 10' 10' 10 
Photon Energy (MeV) 



10" 10" 



10" 10"' 10"' 10' 10° 10° 
Photon Energy (MeV) 



Figure 6. Comparisons of the calculated spectra with the observations for the SNR G23.3— 0.3 at different ages. Synchrotron (dashed 
line), bremsstrahlung (dotted line), Compton (dash-dotted line) of the primary electrons (thin) and the secondary pairs (thick), 
meson decay due to the p-p interaction (black solid line) with radio, X-ray (Tian et al. 2007), and VHE 7-ray with HESS (Aharonian 
et al. 2006b) and MAGIC (Albert et al. 2006) observational results are indicated in the figure. The other parameters are Moj = 1.5 Mq, 
"ISM = 6cm-3, BisM = 5/iG, Kep = 0.2, a distance of 4kpc, and an age of (a) 15000 yr, (b) 30000 yr, (c) 50000 yr, (d) 80000 yr. 
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Figure 7. The comparison of detailed spectra calculated by the 
model with the observations for the SNR G23.3— 0.3 with an age 
of 29800 yr. The other parameters are M^j = 1.5 Mq, nigM = 
1.5cm~'', BisM = SfiG, Kcp = 0.05, and a distance of 4kpc. 
Others are same as those in Fig. [6] 

We now apply the model in this paper to the old SNR 
G23.3— 0.3, which has been observed in radio, X-ray and 
VHE 7-ray bands. Firstly, using the similar process as that 
for the SNR G8. 7-0.1, we find K^p must be about 0.2 for 
the SNR G23.3— 0.3 with parameters, a distance of 4 kpc. 



Moj — 1.5 Mq, msM = 6cm~^, _Bism — 5/iG, and an age of 
several 10* yr. Fig.|6]shows the resulting spectra for different 
ages. For these parameters, the radiative age is 1.43 x 10* 
yr and then the emission in each panel of Fig. [6] is in the 
radiative phase. The observations in X-ray band can be ex- 
plained as the synchrotron and bremsstrahlung radiation of 
the primary electrons for an age of 15000 yr, however, the 
calculated results are significantly high above the "VHE 7-ray 
observations by HESS and MAGIC (see panel (a) in Fig. IS}. 
When the radiative SNR ages, the cut-off energy of the pri- 
mary electrons deceases quickly and the high-energy protons 
undergo both catastrophic energy losses due to the p-p colli- 
sions and energy-dependent diffusion. As a result, the roll-off 
energy of the synchrotron radiation of the primary electrons 
decreases quickly and the number of the TeV photons from 
the SNR reduces too. We can see from panel (d) in Fig. [6] 
that the radio and "VHE 7-ray observations can be explained 
by the model with an age of 80000 yr, however, the observed 
flux in the X-ray band can not be reproduced. To make the 
model results consistent with the flux of the observations in 
the X-ray band, the age of the SNR should not be much 
bigger than the radiative age of the SNR. Furthermore, a 
smaller nisM should be chosen to make the calculation spec- 
trum in the "VHE 7-ray band consistent with the HESS and 
MAGIC observations. 

Finally, nisM is chosen to be 1.5 cm""'^ with other pa- 
rameters an age of 2.98 x 10* yr, Mcj — 1.5 Mq, Bism = 
8fiG, Kcp = 0.05 to model the multi-band non-thermal 
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emission from the SNR G23.3— 0.3 and the resulting spectra 
are shown in Fig. [T] The radiative age for these parameters 
is 2.97 X 10* yr and the roll-off energy of the synchrotron ra- 
diation of the primary electrons is not reduced mostly with 
an age of 2.98 x 10* yr. It is obvious from Fig.[7]that (1) the 
emission from the primary electrons dominates that from 
the secondary pairs in all bands and the emission from 
radio and X-ray band is mainly from the synchrotron ra- 
diation of the primary electrons; (2) the TeV photons are 
produced mainly via the p-p interaction. 



4 CONCLUSIONS AND DISCUSSIONS 

In this paper, by taking into account the evolution of sec- 
ondary pairs produced via the p-p interaction due to 
high-energy protons collide with the ambient matter in 
an SNR, we developed a time-dependant model of multi- 
waveband non-thermal particle and photon spectra both for 
young and for old shell-type SNRs. The primary electrons 
accelerated directly by the shock wave usually dominate the 
secondary pairs produced from the p-p collisions for a 
young SNR, however, the emission from the secondary 
pairs becomes more and more important and eventually sur- 
passes the radiation from the primary electrons as the SNR 
grows old because the source of the primary electrons is cut 
off and the synchrotron radiation loss is significant for elec- 
trons when the SNR is in the radiative phase, whereas the 
secondary pairs can be produced for a long time due to 
the small cross section for the p-p interaction (see Fig. 
Ksp, which can not be directly deduced from the observa- 
tions now and is usually limited by comparison of the model 
results with the multi-band observations for a given SNR, 
is also a crucial parameter to determine whether the con- 
tribution to the final non-thermal emission of the secondary 
e* pairs dominates that of the primary electrons. For an old 
SNR with small K^p, the non-thermal photons with energies 
< 10 MeV are usually produced mainly via the radiations of 
the secondary pairs (see Fig. [2| . 

The model is applied to two old shell-type SNRs, 
G8. 7-0.1 and G23.3-0.3, whose VHE 7-rays are detected 
in the HESS survey of the inner Galaxy. G8.7— 0.1 is a shell- 
type SNR with an age of ~ 15000 yr. The SNR must be in 
the radiative phase to make the calculated results consis- 
tent with the Chandra X-ray observations. Finally, a set of 
parameters, a distance of 6 kpc, an age of 18000 yr, Moj — 
1.5 Mq, nisM = 6cm~^, Kep = 0.15, and Bism = 6/xG is 
used to model the multi-wavelength non-thermal emission 
from the SNR G8. 7-0.1. It can be concluded that (1) the 
radio emission from the SNR is mainly from the synchrotron 
radiation of the primary electrons, whereas the non-thermal 
photons with energies around 0.1 keV primarily come from 
the synchrotron radiation of the secondary pairs pro- 
duced from the p-p collisions; (2) the observed X-ray pho- 
tons with energies ranging from 2 to 10 keV are produced 
mainly via the bremsstrahlung of the primary electrons; (3) 
the TeV photons with energies < 1 TeV are primarily from 
both the bremsstrahlung of the primary electrons and the 
p-p collisions of the primary protons, however, those with 
higher energies are produced mainly through the p-p inter- 
action (see Fig. [5)1. Another old SNR is G23. 3-0.3, which is 
an asymmetric shell-type SNR and spatially coincide with 



HESS J1834-087. Tian et al. (2007) reported new H I ob- 
servations from the VLA Galactic Plane System and a new 
XMM-Newton observation for HESS J1834-087, and argued 
that G23.3-0.3 is an old SNR with a distance of 4 ±0.2 kpc 
and an age of ~ 10^ yr, moreover, a density of 6 cm^^ was 
estimated around the SNR from the VGPs H I line emission 
associated with the SNR. We find that, using nisM ~ 6 cm""^ 
and an age of ~ 10^ yr with other appropriate parame- 
ters, the model result agrees with the radio and VHE 7-ray 
observations, but the calculated X-ray flux is significantly 
lower than the observed one given by XMM-Newton (see 
Fig. A smaller age and a smaller nisM are needed to 
make the calculation result consistent with the multi-band 
observations. Finally, we use parameters a distance of 4 kpc, 
Afoj = 1.5 Mq, an age of 2.98 x 10* yr, nisM = 1.5 cm~^, 
BisM ~ 8 piG, and Kcp = 0.05 to model the multi-band 
non-thermal emission from the SNR G23.3— 0.3. The mod- 
eling results can be consistent with both the radio, VHE 
7-ray observations and the X-ray fiux with these parame- 
ters. Moreover, the emission from the secondary pairs 
is negligible and the multi-wavelength photons with those 
parameters for the SNR mainly come from the radiation of 
the primary particles, and the TeV photons primarily have 
hadronic origin (see Fig. [T}. 

The model in this paper includes the contribution of the 
secondary pairs produced from the p-p collisions to the 
non-thermal emission from SNRs. The secondary pairs 
can be produced for a long time and experience energy losses 
and diffusion, and the contribution of the secondary pairs 
to the final non-thermal emission from an SNR can become 
more and more important and eventually surpasses that of 
the primary electrons when the SNR ages in the radiative 
phase. Whether the contribution of the secondary pairs 
to the final non-thermal emission is prominent or not de- 
pends heavily on the SNR age, A'op and nisM- These pa- 
rameters can be limited by the comparison of the model re- 
sult with the multi-wavelength observations for a given SNR, 
however, the shortage of precise observations and uncertain- 
ties such as the initial eject mass and the initial velocity of 
the shock lessen the reliability of the calculated results. On 
the other hand, for an old SNR with observed TeV emis- 
sion, if the TeV photons primarily have hadronic origin, the 
emission from the secondary pairs is usually prominent. 
The morphology of the resulting spectrum for an SNR in 
the energy range between 10 MeV to 1 TeV directly relates 
to the primary origin of the VHE photons. Fortunately, the 
Gamma Ray Large Area Space Telescope (GLAST) have a 
good sensitivity and resolution in the energy range between 
30 MeV and 300 GeV, and the observations with GLAST 
will give a direct evidence on whether the TeV photons from 
an SNR mainly have hadronic origin or leptonic origin. Ob- 
viously, GLAST will give more limits on the parameters in 
the modeling process with the model in this paper for an 
SNR, and then we can get more insights on the acceleration 
and photon production processes involved in SNRs. 
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